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Antiferromagnetic insulators can become active spintronics components by controlling and de¬ 
tecting their dynamics via spin currents in adjacent metals. This cross-talk occurs via spin-transfer 
and spin-pumping, phenomena that have been predicted to be as strong in antiferromagnets as in 
ferromagnets. Here, we demonstrate that a temperature gradient drives a significant heat flow from 
magnons in antiferromagnetic insulators to electrons in adjacent normal metals. The same coeffi¬ 
cients as in the spin-transfer and spin-pumping processes also determine the thermal conductance. 
However, in contrast to ferromagnets, the heat is not transferred via a spin Seebeck effect which is 
absent in antiferromagnetic insulator-normal metal systems. Instead, the heat is transferred via a 
large staggered spin Seebeck effect. 


PACS numbers: 72.25.Mk,72.20.Pa,73.50.Lw,72.10.Di 

In spintronics, the properties which make antiferro¬ 
magnets markedly different from ferromagnets also make 
them attractive in a more dynamic role. Antiferromag¬ 
nets operate at much higher frequencies and may em¬ 
power Terahertz circuits. They also have no magnetic 
stray fields, which therefore enables denser spintronics 
circuits. For these reasons, antiferromagnets are usually 
passive spintronics components. However, they can play 
a role as active components despite their lack of a macro¬ 
scopic magnetic moment [THIS] and even when they are 
insulating [TUI[T^IT^ . 

We demonstrate that the thermal coupling between 
antiferromagnetic insulators (AFIs) and normal metals 
is relatively strong. The strong thermal coupling facil¬ 
itates several outcomes, can lead to efficient cooling of 
antiferromagnetic spintronics devices, might function as 
heat sensors and can reveal valuable information about 
the high-frequency spin excitations in DC measurements 
that are complicated to extract with other techniques. 

Antiferromagnets can produce pure spin currents as 
large as those produced by ferromagnets. We recently 
showed that spin pumping may be as operative from anti¬ 
ferromagnets as from ferromagnets [13j , in apparent con¬ 
traction to naive intuition. Furthermore, the efficiency 
of spin pumping from antiferromagnets to normal met¬ 
als implies, via Onsager reciprocity relations, that there 
is a considerable spin-transfer torque on antiferromag¬ 
nets from a spin accumulation in adjacent normal metals. 
However, in the absence of external magnetic fields, the 
spin Seebeck effect vanishes [2] . This fact seems to indi¬ 
cate that spins in antiferromagnets decouple from, or are 
only weakly connected to, heat currents and temperature 
gradients in adjacent normal metals. 

To the contrary, we find that the thermal coupling con¬ 
stant is orders of magnitude stronger than its ferromag¬ 
netic counterpart. This radical difference is caused by the 
large exchange field in antiferromagnets that governs the 
heat transfer rather than the much smaller anisotropy 
fields or external magnetic fields in ferromagnets. The 
thermal coupling between antiferromagnetic insulators 


and normal metals is associated with a staggered spin 
Seebeck effect rather than via the spin Seebeck effect. 

Spin caloritronics determines how spins are coupled to 
currents and temperature gradients m- Measurements 
of important thermoelectric properties in ferromagnetic 
insulators, such as the spin Seebeck effect m , are cen¬ 
tral to this field. In the spin Seebeck effect, a tempera¬ 
ture gradient transfers a magnon spin current in a ferro- 
magnet into an itinerant spin current in a normal metal 
m I18j . This process is active even in insulating fer¬ 
romagnets m- The spin Peltier effect is reciprocal to 
the spin Seebeck effect; a heat current generates a spin 
accumulation [iniiiii. These fascinating thermoelectric 
properties can be useful to control the heat flow in spin¬ 
tronics devices and in devices that recycle waste heat. 

In explaining our calculations, we interpret the theories 
on the spin Seebeck effect [HHIHlIua as a combination 
of three mechanisms. First, a precessing magnetization 
can pump a spin current across a ferromagnet-normal 
metal junction |23H25j . Spin-pumping gives rise to an in¬ 
creased magnetization dissipation rate [23l|26l|27|. Sec¬ 
ond, the enhanced dissipation implies that there is also 
an enhanced spin current noise in terms of a fluctuating 
spin-transfer torque |29j . In equilibrium, there is no ther¬ 
mal bias and the DC spin current vanishes because the 
temperature-driven spin pumping and a fluctuating spin- 
transfer torque exactly compensate each other. Third, a 
temperature difference alters this balance and causes a 
net spin current [muHiiiui. 

In this picture, to compute the heat transfer between 
AFIs and normal metals, we first establish the fiuctu- 
ating spin-transfer and staggered spin-transfer in such 
hybrid systems. Both quantum and thermal fluctuations 
are required to determine the magnon occupations. Sub¬ 
sequently, we use these results to define the thermal gra¬ 
dient driven (staggered) spin currents, which we then use 
to evaluate the rate of change of the heat in the AFI. We 
focus on insulating antiferromagnets where the transport 
properties are magnon-driven. Generalizations to con¬ 
ducting antiferromagnets are straightforward. 
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Figure 1: An antiferromagnet insulator (AFI) sandwiched be¬ 
tween two normal metals A^i and N 2 - The left normal metal 
(Ail) is a good spin sink. The right normal metal {N 2 ) has no 
spin memory loss. A heat current between the normal metals 
flows in response to an applied temperature gradient across 
the AFL The cross section is A and d is the AFI thickness. 
The heat flow Iq is along the longitudinal coordinate x. 


the spin sink and the subscript v denotes the product of 
either the sub-index n or m and a Cartesian component 
X, y, or z. The bulk torques arise from the magnon- 
phonon interaction. At the AFI-N interface, the torques 
are governed by spin-pumping induced by the coupling of 
the magnetic moments to itinerant electrons in adjacent 
normal metals. In finding the torques, we introduce dissi¬ 
pation in a similar way as in Ref. [7] and further generalize 
this description to include quantum and thermal fluctu¬ 
ations. The resulting fluctuation-dissipation torques are 

T„ = (hm — am) X n -I- (h„ — an) x m, (4a) 

Tm = (h„ — an) X n -b (h^ — am) x m, (4b) 


We model the AFI as a two sub-lattice system with 
spatio-temporal magnetizations Mi and M 2 . The dy¬ 
namics are described by the staggered magnetizations 
L = Ml — M 2 = Ln and the magnetization M = 
Ml -b M 2 = Lm. These fields satisfy the constraints 
n^ -b m^ = 1 and n • m = 0. In equilibrium, the stag¬ 
gered field is homogeneous and constant in time, |L| = L, 
and the magnetization vanishes, M = 0, i.e. n^ = 1 and 
m = 0. We consider an easy-axis AFI that is described 
by the free energy F = J dr [/(r) + /s(r)], where the en¬ 
ergy density is 


/=- 

7 


- n^) - + nl) 


( 1 ) 


with the exchange frequency uje and the smaller 
anisotropy frequency wa, ^ When n and m 

spatially vary, the stiffness contributions are 

Is = X! -b {Xmdiinf] , (2) 

i—x,y,z 

where and are exchange lengths associated with 
n and m, respectively. The dynamic equations are: 


n = cJrrj X n + u;r7 X m - 


m. = (jJri X n -h X HI - 


(3a) 

(3b) 


where the effective fields are a)„ = —('y/L)Sf/Sn and 
ujjn = —{'y/L)6f/Sm. In Eqs. (3a) and (3b), the dissipa¬ 
tion and fluctuation torques Tm and are essential to 
describe spin caloritronics effects. 

We consider a thin-film AFI of thickness d sandwiched 
between two normal metals, the left one of which is a 
good spin sink (e.g. Pt), and the right one has little or 
no spin memory loss (e.g. Cu), see Fig. We assume 
planar AFI-normal metal interfaces of cross section A. 
The coordinate r = (x, p) is decomposed into a perpen¬ 
dicular coordinate x (0 < x < d) and the 2D in-plane 
coordinates p. The fluctuation-dissipation torques have 
bulk and (spin-pumping induced) interface contributions, 
T,y(r) = Ti^\r) + 5{x — xi)Tlf\p), where xj = 0+ is in¬ 
finitesimally near the interface on the AFI side close to 


for both bulk and interface contributions where we have 
suppressed the superscript [(&) or (p)] in t, h, and a. 
The bulk Gilbert damping constant is a*^^^ is a 

measure of the spin-pumping-induced enhanced dissipa¬ 
tion; for homogenous macrospin excitations the enhanced 
damping constant is a^^^ + /d [25]. 

The fluctuation-dissipation theorem implies the exis¬ 
tence of the fluctuating forces and h„. The average 
of the independent fluctuating forces and h„ as well 
as all their fluctuations vanish, except 

(hy(pt)/zW(pr)) = , (sa) 

(hW(rt)hW(r't')) = . (5b) 


The correlation function R{t, T) depends on the (local) 
temperature. As demonstrated for ferromagnets in Ref. 
1291 the spin-current fluctuations associated with spin¬ 
pumping depend on the temperature in the normal metal 
close to the interface, Ti. We posit that the one-to-one 
correspondence between spin-pumping in ferromagnets 
and antiferromagnets m implies that the spin-current 
fluctuations in antiferromagnets obeys the same relation¬ 
ship, as in Eq. (5aI. In the bulk of the AFI, the phonon- 
induced fluctuations associated with the bulk Gilbert 
damping depends on the temperature profile in the anti¬ 
ferromagnet Ta{x), as in Eq. (5b). The correlation func¬ 
tion only describes white noise in the (classical) high- 
temperature limit, R{t,T) « 27rkBT6{t). However, for 
our purposes, the magnon occupations also depend on 
the low-temperature quantum limit of the fluctuations, 
which we describe after Eq. (16). 

The effective fields determined by Eqs. ([^ and ([^ are 


a;„ = -b a;A(n • z)z + a;A(A„V)^n, (6a) 

=-W£;m-bWA(m-i)z-bWA(AmV)^m. (6b) 

We focus on the rate of change of the heat {Iq) from 
the AEI into the normal metal. Evaluating the rate of 
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change of the free energy due to the (spin-pumping) inter¬ 
face fluctuation-dissipation torques, Iq = f dplujn ■ 
-|- iJrn ■ and using the continuity equations for 
the spin and staggered spin currents below in Eqs. 
we find 


LAd f dp 




dn_ 
' ddx 


( 7 ) 


where we use a circular basis so that n± = ± iUy and 


m± = rrix ± irriy are first order corrections with respect 
to the equilibrium configuration n = z and m = 0. Iq 
has units of energy per time. 


Next, we Fourier transform in the transverse coor¬ 
dinate p and time t so that any function c(a;, p, t) = 
/ c?a;c(a;, q, w) expi(a;t — q • p). Using Eq. (j^, the 
linearized dynamic equations of motion become 


-t-w^ [l -I- A^(q^ - 5^)] -I- 2uje'^ m± = ±ojn± + , 

-f UJA [(1 + A^(q^ - dl)] | h± = ±a;m± -h h^ni{x). 


( 8 a) 

( 8 b) 


In the coupled dynamic equations Q, the stiffness contri¬ 
butions can be interpreted as arising from the conti¬ 
nuity equations for the staggered field and the magnetiza¬ 
tions, (9tn)s = 0 and (9tm)s = 0 . 

In linear response, the staggered spin current and spin 
currents along the ^-direction are jn,x = ^ 

and im,x = ^ 9a;n. The boundary conditions for 

the linearized equation of motion are obtained by in¬ 
tegrating the dynamic equations ra across the AFI-N 
interface. This results in the continuity of the spin and 
staggered spin currents in linearized forms at a; = 0 : 


dm± . ( ) 


= iu!a''^^rh± — , 


dx 


waA: 


2 9h± . (p) 


dx 


= iuja'^^^h± — h^n'!± ■ 


(9a) 

(9b) 


Similarly, at a; = d, there is no loss of currents and 
the boundary conditions are CLi^A^9m±/9x = 0 and 
WA^'iidfi±ldx = 0 . 

In typical antiferromagnets, uje is much larger than 
all other energy scales and we may employ the so-called 
exchange approximation. This implies that we may dis¬ 


regard smaller terms in the equation of motion ( 8 a) so 


that it greatly simplifies to 7 ti± = Lijh±l2u}E- By insert¬ 
ing this relation into Eq. ( 8 b), we find the equation of 


motion in the exchange approximation 


n± 


K{qi + di)n± = -^ 


UJA 


which can be solved with the boundary conditions of Eq. 


(9b). In the exchange approximation, to the lowest order 


in the dissipation, we have introduced the longitudinal 
wavenumber q^- The complex wavenumber q^ is implic¬ 
itly defined via the relation w = ujr + where the 

bulk resonance frequency and the bulk lifetime are deter¬ 


mined by 

ujj^ = 2 uJaUJe{^ + ^niQx + dy + 9z)) ; (H^) 

1/tW ws. (lib) 

The central results we will obtain can be interpreted 
in terms of the eigenstates with the associated eigen- 
frequencies and lifetimes in a thin film antiferromag- 
net. The eigenstates are determined by expressing fi± = 
A±e'xp{iqxx) + B± exp {—iq^x) in Eq. (10) when the 
right hand side (the fluctuations) vanishes. The only 
non-trivial solution that satisfies both the boundary con¬ 


ditions of Eqs. (9b) at a; = 0 (with no fluctuations) and 


dh±/dx = 0 at a; = d is determined by the secular equa¬ 
tion s[qx) = 0 , where 


5(dx) = 


qxK^^A 


tan (qxd) — i- 


( 12 ) 


doj d 

In the absence of spin-pumping and bulk damping, the 
solutions of s{qx) = 0 are standing waves where qx = 
Nir/d and N is an integral number. When spin-pumping 
is weak, the second term in Eq. (12) is small and the 


solutions of s{qx) = 0 can be expanded around the solu¬ 
tions obtained in the absence of spin-pumping. For the 
higher modes, when iV 7 ^ 0, we expand the wavevector qx 
to the first order in the deviations from Nn/d and insert 
the resulting imaginary part of the wavevector into the 


lifetime For = 0, we carry out a second-order 

( 10 ) 

expansion in terms of the small parameter around 0 


dispersion relation of Eq. (11a) to find the spin-pumping 


and insert this result into the dispersion of Eq. (11a) to 

(p) ' -' 

find the lifetime tg . We compute that 


i/t?’ = 


1 


= 2 



(13a) 

-^LJE ■ 
a 

(13b) 
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In a striking contrast to ferromagnets, the spin-pumping 
induced scattering rate is proportional to 

the exchange energy rather than the spin-wave eigenfre- 
quency. This is in agreement with the behavior for the 
bulk lifetime (lib), and, as we will see, implies a large 
heat conductance between APIs and normal metals. In¬ 
terestingly, we hnd that the spin-pumping induced relax¬ 
ation rate of the higher modes is twice as large as the uni¬ 
form, but independent of the transverse (2D) wavevector 
q. This ratio agrees with our previous result for the spin¬ 
pumping induced ratio in thin film ferromagnets and can 
be used to distinguish the spin-wave modes [ 50 ] . 

Next, we solve the linearized dynamic equation ([^ 
with the fluctuating bulk forces and subject to the bound¬ 
ary condition (9b) where the fluctuating spin-pumping 
induced forces appear. To compute the rate of change of 
the heat, we represent the solution at x = 0 as ffi^ = 


+ Jo dxx^J^l(x)hl^l(x) and ujA^ldxn+jd = 




(p) lip) 

Xn^+hnl + fo dxx'in'+ix)h^^l{x). We hnd that = 
-l/{2duJEsiqx)), Xm+ = Xm+ cos qa;id - x)/cos qa;d, 
= -qj;X'^uJAtanq^d/d‘^ujs{qj;), and XnL = 


X. 


(p) 
n' + 


cos qa:{d — x)/d^s{qx) cos qxd. 

We evaluate the variance of the huctuating forces 
and hnd the rate of change of the heat, Iq = 
(2dw_E/7r)ImX)q/.!^ -|- 77^^'^], where the spin¬ 

pumping and bulk contributions are 


Vq'’ = xl^'’+iXn^+)*a^^^Ri‘^,Ti ), (14a) 

= f dxx^^+[x)[XnK^x)ya^'^'>R{uj,TA{x)). (14b) 
Jo 


In equilibrium, Ta{x) = Ti, the rate of change of the heat 
vanishes, Iq = 0, as expected. In linear response, the 
temperature varies linearly in the AFl so that Ta{x) = 
Ti -\- (T2 — Ti)x/d. We then compute that the rate of 
change of the heat is Iq = (T2 — T'i)kq, where 


KQ 


= E 




(p) 


duj 


V ,dR{u:,T) 


Qy,Qz 


‘2TT\s{qx)\' 


yilx) 


dT 


(15) 


and Cilx) = ‘2 Jo ^. By following the 

same methods, we also compute that the temperature- 
driven spin current, i.e. the spin Seebeck effect, van¬ 
ishes, in agreement with Ref. 1141 However, we find that 
the temperature-driven staggered spin current is finite. 
The heat is therefore transferred between the AFl and 
the normal metal by means of a staggered spin current. 

By comparing the equilibrium expectation value of the 
spin-wave free energy with the quantum-mechanical re¬ 
sult for a magnon gas or, alternatively, by using the 
fluctuation-dissipation theorem represented by Eq. (4.9) 
in Ref. [STJ we identify that the correlations function 
R{uj, T) represents the mean energy at the temperature 
T of an oscillator at natural frequency w, R{uj, T) = 
+ hLuf{uj,T), where f{uj,T) is the Bose-Einstein dis¬ 
tribution function. At high temperatures we thus ob¬ 
tain the classical R « ksT, as expected. However, the 
complete frequency dependence of R{ijJ,T) is essential 
in determining the thermally active magnon modes and 
quantifying the heat conductance kq. 


In the limit of low bulk damping and spin-pumping induced damping, we can expand the poles of the denominator 
of Eq. (15) around the spin-wave resonance qxd = mr in a similar way as in Ref. 1221 This results in an intuitive 
expression: 


1 

dwDw(a;)fia;{/(a;,Ti)[l-/(a;,r2)]-/(a;,r2)[l-/(a;,ri)]} . (16) 

N=0 


The heat current that flows between the nor¬ 
mal metals via the antiferromagnet, at each fre¬ 
quency, is proportional to the spin-pumping in¬ 
duced spin-wave relaxation rate l/t^\ the mode- 
dependent density of states, Dn{uj) = J2qy ~ 


{2ujaoje [1 + + 92)^] }^^^)- Furthermore, 

the heat current is determined by the Bose-Einstein oc¬ 
cupation of the magnons and the electron-hole pairs in 
the normal metal. This expression (16) reveals that the 
thermal coupling between normal metals and AFIs is 
relatively strong. The heat current is proportional to 
the spin-pumping induced spin-wave scattering rates that 


are proportional to the exchange energy and the Gilbert 
damping coefficient. At high temperature, we find h = 
ATrykBTiykB{T2-Tya(-Py{15V2AiyyU^h^/y,where 
Aex = fujOA^J is the exchange stiffness. 

Phonons also mediate heat currents between AFIs and 
normal metals. Experimentally, the magnon-induced 
heat current we predict here can be separated from the 
phonon heat current by the different temperature and 
length dependence. Eor instance, different measurements 
of normal metals with strong and weak spin-orbit cou¬ 
pling can be compared. Additionally, one can use an 
external magnetic held to change the magnon dispersion 
and consequently the spin-wave density of states govern- 
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ing magnon-induced heat current of Eq. (16). 

In conclusion, we demonstrated a strong thermal cou¬ 
pling between antiferromagnetic insulators and normal 
metals. The heat is carried in terms of a staggered spin 
current. 
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Norway, project number 216700. 
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